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A comparison of calculated and experimental Qa values of superheavy even-even nuclei and a few 
selected odd-A'^ nuclei is presented in the framework of the relativistic mean-field model with the 
parameterization NL-Z2. Blocking effects are found to be important for a proper description of Qa 
of odd mass nuclei. The model gives a good overall description of the available experimental data. 
The mass and charge assignment of the recently measured decay chains from Dubna and Berkeley is 
in agreement with the predictions of the model. The analysis of the new data does not allow a final 
conclusion about the location of the expected island of spherical doubly-magic superheavy nuclei. 
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In the last few years, the synthesis of superheavy nuclei 
with Z = 110-112 at GSI (Darmstadt) and JINR Dubna 
- see Q and references therein - has renewed the interest 
in the properties of superheavy nuclei. These are by def- 
inition the nuclei with Z > 100 which have a negligible 
fission barrier and are stabilized by shell effects only. The 
ultimate goal is to reach an expected 'island of stability' 
located around the next spherical doubly-magic nucleus 
which was predicted to be 184II4 thirty years ago [^]. 

Recently the discovery of new rather neutron-rich iso- 
topes of the elements Z = 112 |Q] and Z = 114 ^ was 
reported from JINR Dubna, while an experiment at 
Berkeley led to the synthesis of even heavier nuclei with 
Z = 118 Q. While earlier superheavy nuclei could be un- 
ambiguously identified by their a-decay chains leading to 
already known nuclei, the decay chains of the new-found 
superheavy nuclei cannot be linked to any nuclei known 
so far. Their identification relies on comparison with the- 
oretical models. 

The experimental progress is accompanied by a large- 
scale investigation of superheavy nuclei with the latest 
nuclear mean-field models. While refined macroscopic- 
microscopic (MM) models like the Yukawa-plus- 
exponential model with Woods-Saxon single-particle 
potentials (YPE-t-WS) or the finite-range droplet 
model with Folded- Yukawa single-particle potentials 
(FRDM-HFY) confirm the older prediction of UllU 
for the next spherical doubly-magic nucleus, most self- 
consistent models shift that property to the more proton- 
rich mi20 or even 111126 |o|-^. The conflicting pre- 
dictions have several reasons. MM models have very re- 
stricted degrees of freedom of the radial density distri- 
bution and the shape of the single-particle potentials, 
which becomes visible in superheavy nuclei where this 
hinders the occurrence of proton shells at Z = 120 or 
Z = 126 |l|,|l|,|l|l. Another source for different extrapo- 
lation to superheavy nuclei among the models is the spin- 
orbit interaction. While the self-consistent relativistic 



mean-field model (RMF) naturally incorporates the nu- 
clear spin-orbit interaction (which is a purely relativistic 
effect), it has to be put in by hand into all non- relativistic 
models, self-consistent ones using Skyrme (SHF) interac- 
tions and MM models. Surprisingly all non-relativistic 
models - which have one or several additional parame- 
ters to explicitly adjust the experimental spin-orbit split- 
tings - perform not as good in this respect as the RMF 
which has no parameters adjusted to single-particle data 
at all. The non-relativistic models overestimate spin- 
orbit splittings with increasing mass number which might 
have some impact on their actual predictions for shell 
closures in the superheavy region |l4| . Superheavy nu- 
clei with their large density of single-particle states are a 
sensitive probe for models of nuclear shell structure and 
can be used to discriminate among models which describe 
stable nuclei with comparable quality. 

Macroscopic-microscopic models - which have been the 
favorite tool to describe superheavy nuclei for a long time 
- provide a very good description of masses through- 
out the chart of nuclei, but they rely on preconceived 
knowledge about the densities and single-particle poten- 
tials which fades away when going towards the limits of 
stability. Although self-consistent models have not yet 
reached the overall performance of MM models, they in- 
corporate polarization effects on the density distribution 
which might be crucial for a proper description of super- 
heavy nuclei where the strong Coulomb field pushes the 
protons to the nuclear surface pUJT^ . 

Recently Cwiok et al. have reported a detailed compar- 
ison of the new experimental data with SHF calculations 
employing the interaction SLy4 |16|. Although the cal- 
culated Qa values show some deviations from the mea- 
sured values, the analysis confirms the assignment of the 
mass and charge number of the new nuclides. SLy4 shifts 
the island of stability towards very high charge numbers 
around 134 126 |11-14|. It is the aim of this contribution 
to compare the experimental data with predictions of the 
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FIG. 1. Qa values of even-even nuclei in the superheavy 
region calculated with NL-Z2 (open circles) compared with 
experimental values for even-even nuclei (filled diamonds) 
where available and odd-Af isotopes (open diamonds) of the 
heaviest even Z elements. Solid lines connect nuclei with the 
same proton number Z, dotted fines nuclei in a-decay chains. 
Including data for nuclei with Z > 116 leads to overlapping 
curves, therefore those are omitted in the plot. 

RMF model obtained with the recent parameteriza- 
tion NL-Z2 by P.-G. Reinhard This particular force 
provides the best overall description of binding energies of 
superheavy nuclei among the current parameterizations 
of the RMF ||l|,|l|| and with ?7il20 gives an alternative 
prediction for the next doubly magic nucleus which is 
quite close to the upper end of the new Z = 118 chain. 

The calculations are performed on an axially symmet- 
ric grid assuming reflection symmetry. Nuclei with odd 
mass number are calculated in a self-consistent block- 
ing approximation taking into account time-odd contri- 
butions to the Dirac equation as described in Pair- 
ing correlations are treated within the BCS scheme using 
a delta pairing force, see for details. There remains 
a small uncertainty of the calculated Qa values on the 
order of a few 100 keV from a possible variation of pair- 
ing recipes (i.e. particle-number projection and choice of 
pairing interaction), especially in odd nuclei. The pairing 
strength is fitted to pairing gaps calculated in spherical 
blocking approximation. Taking polarization effects in 
odd nuclei into account during the fit gives larger pairing 
strengths |2^j2l[], but this affects mainly the odd-even 
staggering, while the Qa of heavy nuclei are rather ro- 
bust. The correction for spurious center-of-mass motion 
is performed a posteriori as done in the original fit of 
NL-Z2, while corrections for spurious rotational and vi- 
brational motions are neglected (as it is done in all other 
studies of superheavy nuclei) which might affect Qa val- 
ues in some cases to the order of 1 MeV and tend to 
diminish shell effects p^ . 

A first test of the reliability of the RMF to describe 
the new data is to check its performance for known Qa 
values of even-even nuclei, see Figure 0. Owing to the 
lack of data for even-even nuclei the Qa of odd-iV nu- 



clei are drawn beyond Z = 104. They have to be han- 
dled very carefully, some of these Qa values might corre- 
spond to transitions involving excited states, and due to 
the blocking effects the ground-state-to-ground state Qa 
values might deviate on the order of 500 keV from the 
systematics of the Qa values of even-even nuclei. NL- 
Z2 gives a good overall description of the data except 
for some nuclei around Rfio4 where NL-Z2 overestimates 
a deformed proton sub-shell closure while the deformed 
N = 152 shell is underestimated. The latter is a prob- 
lem from which virtually all self-consistent models suffer 
1 13,^. The particular Skyrme interaction SLy4 used in 
1 16 1 performs better in the region around Z = 104, but 
it is to be noted that this interaction has small errors in 
the isotopic and isotonic mass systematics of these nuclei 
||l3|,|8|, which cancel when calculating Qa values. 

To get an impression of the shell structure (as predicted 
by NL-Z2) of nuclei in the new-found region. Fig. |^ shows 
the single-particle spectra of the even-even nuclei with 
one neutron less than those in the decay chain of 118. 
The nuclei close to the Z = 120 shell closure come out 
spherical including ^y|ll8, while all lighter nuclei have 
deformed ground states. In addition to the spherical 
Z = 120 shell, the proton spectra show several deformed 
shell closures indicated in the plot. In contrast to MM 
models and some SHF interactions, however, the proton 
spectra show no preferred (deformed) shell closure for the 
light nuclei in this chain, but a region of small level den- 
sity between Z = 104 and Z = 110. The single-particle 
spectra of the neutrons are much denser with a spheri- 
cal iV = 172 shell and deformed shell closures at iV = 174 
and N — 162. The latter ones are also predicted by other 
models [|^,@,16| and N = 162 is already confirmed exper- 
imentally |^[~^The single-particle spectra of nuclei close 
to the decay chain of 175 114 look similar except for the 
fact that all nuclei there are deformed. 

Three a decay chains of odd-mass nuclei are discussed 
in the following: the one starting with ig5ll2 measured 
at GSI in 1994 ||2^ serving as a testing ground for the 
performance of NL-Z2 when describing odd-mass number 
nuclei, and the new data assigned to f7|ll4 and 275II8. 
The experimental data were obtained from two, one and 
three events respectively. Fig ^ compares calculated and 
measured Qa values of these selected chains. Values cal- 
culated with the FRDM+FY model 0, the YPE-fWS 
model and the SIIF-|-SLy4 model JM] are given for 
comparison. As can be seen from Fig. |pthe quasiparti- 
cle spectra of odd-iV nuclei in this region are very dense. 
Each single-quasiparticle state is the band head of a ro- 
tational band with transition energies from the first ex- 
ited state to the band head below 50 keV which further 
increases the level density significantly. It cannot be ex- 
pected that the synthesis of these nuclei and their decay 
chains lead always to ground states but very often to ex- 
cited states. Experimental evidence for that was found in 
the decay of fg|ll2, each of the two nuclei identified so far 
p3[ has decayed through a different state in fggllO. The 
branch with large Qa in ill HO presumably connects low- 
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lying states in both nuclei while the branch with the low 
Qa goes through a highly excited state in that nucleus. 
This demonstrates a fundamental problem when com- 
paring measured and calculated Qa of decay chains with 
a few identified events only: the low statistics does not 
allow to identify the whole spectrum of possible a transi- 
tions. Detailed information on the a-dccay fine structure 
of superheavy nuclei is available for a few nuclides up to 
Z = 104 only |^5|. Some guidance when transitions to 
or between excited states might be favored is given by 
the fact that among transitions with similar values 
those between states with the same quantum numbers 
are favored. 

In view of the remaining uncertainties NL-Z2 gives a 




Neutron Number N 

FIG. 2. Single-particle spectrum of the protons (upper 
panel), neutrons (middle panel) and the relative quadrupole 
and hexadecapole moments given by (3i — 'iivir'^Yio) / {iAR^) 
with R = 1.2 A^^^ fm (lower panel) for the a-decay chain 
IfeUO ??Ill8 ■■■ ^ ?6oR'fi04 of even-even nuclei adja- 
cent to the a-decay chain of 175II8. For the deformed nuclei 
at small neutron number the states are labeled with angu- 
lar-momentum projection and parity, while for the spherical 
nuclei at large neutron number the total angular momentum 
and parity are given. The full (dotted) lines in the upper 
and middle panels denote positive (negative) parity states, 
the dashed lines the Fermi energies. 
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FIG. 3. Comparison of experimental and calculated Qa 
values for the decay chains of lllll2, ?7ill4, and frillS, in 
the latter two cases following the mass and charge assignment 
of the experimental groups. In the 165II2 chain two distinct 
branches leading through different states of the intermediate 
nuclei are known. The calculated values from RMF-I-NL-Z2 
and SHF-|-SLy4 connect the lowest states with positive parity 
in all cases (in the new chains only 175 114 and igyllO are 
predicted by NL-Z2 to have ground states with negative parity 
in agreement with @), while the FRDM-fFY and YPE-(-WS 
data are ground state to ground state values. 

very good description of the heavy nuclei in the decay 
chain of 155 112 above N > 157 and reproduces the large 
jump in Qa between lllnO and IfilOS caused by the 
deformed iV = 162 shell closure, which cannot be seen 
in the predictions of the FRDM-I-FY model. For small 
neutron numbers the Qa calculated with NL-Z2 are too 
small as in the case of the even-even nuclei discussed 
above. NL-Z2 gives also a rather good description of the 
a-decay chain of ^^1114. The missing kink in the calcu- 
lated values can be explained assuming that the decay 
of ?75ll4 leads to one of the numerous low- lying excited 
states in ^73112. 

The decay chain of f 7^118 leads through a region with 
seemingly constant total shell effects, while most mean- 
field models predict several spherical or deformed shell 
closures in this region. The FRDM-I-FY data show a 
pronounced sheU effect for 17^114, while the SHF-|-SLy4 
results show one for 173 116 which is related to deformed 
shells at Z = 116 and TV = 174 ||l6|. NL-Z2 predicts a 
deformed = 174 shell as well but a deformed Z — 114 
proton shell (see Fig. ||) that leads to a broad plateau 
at smaller atomic numbers. The predictions of the 
YPE-I-WS model of 0,|2|l follow qualitatively the exper- 
imental data for the heavy nuclei in this chain but show 
a shell effect for ^57 110. In contrast to the other models 
discussed here the YPE+WS model predicts the nuclei at 
the upper end of this decay chain to be oblate deformed 

0- 

The mass and charge assignment of the Dubna and 
Berkeley data are based on theoretical models. While 
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FIG. 4. Comparison of experimental and calculated Qa 
values for the a decay chain of the Z = 119, nuclei recently 
measured in Berlieley assuming other neutron numbers of the 
initial nuclei. The Qa values connect the lowest states with 
positive parity. 



N = 112 predicted by the RMF are relatively weak and 
do not significantly change the systematics of Qa values 
of odd-A^ nuclei at the upper end of the Z = 118 chain 
measured in Berkeley. 
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North Carohna), and DE-AC05-96OR22464 with Lock- 
heed Martin Energy Research Corp. (Oak Ridge National 
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the proton number is rather certain in both cases, the 
neutron number A^ might also be smaller. Fig ^ shows 
the Qa values measured in Berkeley in comparison with 
calculated Qa values for decay chains with varied neu- 
tron number of the initial nucleus. The heavy nuclei 
in all decay chains have very similar Qa values, while 
there are visible differences for the lower end of the de- 
cay chain. It can be clearly seen that Qa of odd nuclei 
do not necessarily follow exactly the trend of the Qa of 
the adjacent even-even nuclei. Since the 2n and higher 
reaction channels are energetically forbidden and statis- 
tically suppressed, an interpretation of the Berkeley data 
in terms of the a-decay chain of 293 118 is in agreement 
with predictions of the RMF with NL-Z2. An interest- 
ing feature of the NL-Z2 results is that Qa values cal- 
culated between the lowest negative-parity quasiparticle 
states starting with 29^ 118 or 29III8 follow very closely 
the experimental data, but several of the intermediate 
nuclei can be expected to 7 decay to states in rotational 
bands build on lower-lying quasiparticle states with pos- 
itive parity which prevents a long a-decay chain between 
negative-parity states (see also ) . 

To summarize: The RMF with the parameterization 
NL-Z2 gives a reasonable description of the Qa values 
of known superheavy nuclei with an overall quality com- 
parable to other mean-field models that predict differ- 
ent spherical magic numbers. On the basis of the avail- 
able experimental data, a decision about the location of 
the 'island of stability' cannot be made. The analysis of 
the new data from Dubna and Berkeley in terms of the 
RMF model suggests that the relatively small Qa val- 
ues and corresponding long half-lives are caused by de- 
formed Z = 114 and A^ = 174 shell closures rather than 
the vicinity of the (potential) doubly magic 1I4II4 pre- 
dicted by some other models, but a spherical Z = llA 
shell, restricted to higher neutron numbers, cannot be 
excluded. The spherical shell closures a.t Z = 120 and 
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